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(XRFs) poses problems both for a power-law universal jet model where the energy per solid angle 
decreases as the inverse square of the angle respect to the jet axis, and for a conical jet model with a 
uniform energy density within the jet beam and a sharp energy cut-off at the jet edge. Here we show 
that the current GRB-XRF prompt emission/afterglow data can be understood in terms of a picture in 
which the GRB-XRF jets are quasi-universal and structured, with a Gaussian-like or similar structure, 
i.e., one where the jet has a characteristic angle, with a mild variation of energy inside and a rapid (e.g. 
exponential) decrease of energy outside of it. A Monte Carlo simulation shows that the current data 
is compatible with such a quasi-universal Gaussian jet with a typical opening angle of 5.7^2^ degrees, 
and with a standard jet energy of about log{Ej/l erg) ~ 51.1 ± 0.3. According to this model, the true- 
to-observed number ratio of the whole GRB-XRF population is about 14 with the current instrumental 
sensitivity. 

Subject headings: Gamma ray bursts - X-rays - gamma rays 



1. INTRODUCTION 

The achromatic steepening in the afterglow lightcurves 
of some GRBs as well as energy budget arguments suggest 
that GRBs are produced by coUimated jets (Rhoads 1999; 
Kulkarni et al. 1999; Harrison et al. 2001). This sugges- 
tion receives indirect support from the intriguing fact that 
a geometry-corrected jet energy appears to be standard 
(Frail et al. 2001; Panaitescu & Kumar 2001; Bloom, Frail 
& Kulkarni 2003; Berger, Kulkarni & Frail 2003a). There 
are two straightforward interpretations of this fact. One is 
that different GRBs collimate the same total energy into 
different angular openings (ranging from 1 to 30 degrees) 
(Rhoads 1999; Frail et al. 2001). Another is that aU GRBs 
have a universal jet shape with a varying energy per solid 
angle which is of the form e{0) oc 6*"^, where 9 is the polar 
angle respect to the jet axis (Rossi, Lazzati & Rees 2002; 
Zhang & Meszaros 2002a). In principle, an individual 
GRB could also have other jet structure, e.g. a Gaussian 
or even arbitrary function (Zhang & Meszaros 2002a). 
When the jet parameters are allowed to have some dis- 
persion around the mean values (which is demanded by 
the data), one has a "quasi- universal" jet structure, and a 
power-law jet structure is no longer a pre-requisite (Lloyd- 
Ronning, Dai & Zhang 2003). 

Recently the so-called X-ray flashes (XRFs, Heise et al. 
2003; Kippen et al. 2003) are identified and proposed as a 
natural extension of GRBs into a softer and fainter regime. 
Another intriguing empirical fact is that the GRB-XRF 
spectral break energy {Ep) in the cosmic rest frame ap- 
pears to be correlated with the "isotropic-equivalent" en- 



ergy (or luminosity) of the explosion (Amati et al. 2002; 
Sakamoto et al. 2003; Lamb, Donaghy & Graziani 2003; 
Lloyd et al. 2000) according to 
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Although the data sample in the XRF regime is currently 
very small (because the faintness of XRFs hinders the 
detection of their possible afterglows and the measure- 
ments of their redshifts), this intriguing result, if confirmed 
by further future data, strongly suggests that GRBs and 
XRFs are related events. A successful model should there- 
fore be able to interpret both GRBs and XRFs within a 
unified framework. It is worth noticing that GRB 980425 
has a very low Eiso (~~ 10'*^ ergs) but a relative large Ep 
(~ 70 keV), which is apparently inconsistent with eq.(l). 
Because it is atypical in other aspects as well, we do not 
interpret it within the current framework. Since there is 
no unified terminology, in this paper we call the combined 
population of GRBs and XRFs as GRB-XRF. 

In this Letter, we do not discuss the nature of the 
Ep oc {Eisoy^'^ empirical law^. Rather we conjecture that 
it is a universal law for the majority of GRB- XRFs, and 
evaluate its implications for various GRB-XRF jet models. 
We will show that a quasi-universal Gaussian-like struc- 
tured jet model is compatible with the current prompt 
emission and afterglow data of GRBs and XRFs. This is 
consistent with earlier suggestions that GRBs would be 
seen as X-ray transients when viewed at large angles (e.g. 
Meszaros & Rees 1997; Woosley, Eastman & Schmidt 1999 
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and subsequent work). 

2. CONSTRAINTS ON OTHER JET MODELS 

A direct implication of the Ep oc {EiaoY^"^ law is that it 
poses important constraints both on the specific "univer- 
sal" structured jot model that assumes a power-law struc- 
ture with index 2 (or k ~ 2 power-law model, Rossi et al. 
2002; Zhang & Meszaros 2002a), and on the conventional 
"uniform" jet model (Rhoads 1999; Frail et al. 2001)^. 

The constraint on the k = 2 power-law universal model 
was raised by Lamb et al. (2003). The argument is as 
follows: In the k = 2 model, we have Eiso oc 6~'^, where 
9y is the viewing angle (since e{9) is equivalent to Eiso in 
the structured jet model). This implies Ep cx 9^^. Since 
the Ep values range from >,300 kcV in GIlBs all the way 
down to <3 keV in XRFs, the viewing angles of XRFs 
need to be two orders of magnitude larger than those of 
the GRBs. Even if XRFs originate in standard jets viewed 
at the equator {9^ ^ 90°), all GRBs should have viewing 
angles less than 1°. The probability for a viewing angle 6y 
is P{6y)d9y (X sm9yd9y, thus a A; = 2 model greatly over- 
predicts the number of XRFs. This is in sharp contrast 
with the observations, which indicate that XRFs repre- 
sent approximately 1/3 of the total GRB-XRF population 
(Lamb et al. 2003). 

The above argument may be regarded as supporting 
a uniform jet model (Lamb et al., 2003), where the ra- 
diation is seen "on-bcam" . However, it docs not ad- 
dress the afterglow data such as light-curve breaks, and 
leads to energetic and progenitor number inconsisten- 
cies. This is because in this model, even if XRFs cor- 
respond to isotropic events, GRBs have to still be jets 
with typical opening angles less than 1°. In the stan- 
dard afterglow model, the bulk Lorentz factor evolves as 
Tit) ~ 6(i;52/n)V8(t/l day)-3/8(l + zf/\ where E52 is 
the isotropic kinetic energy of the fireball, n is the inter- 
stellar medium density, t is the observer's time, and z is 
the redshift. Taking the standard view that the jet break 
time tj (around days) corresponds the epoch of l/r(f) = 6j 
(Rhoads 1999), we get the constraint 
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which requires an extremely large kinetic energy, or an ex- 
tremely low medium density, or both. This is in sharp 
contrast with current afterglow analyses (Freedman & 
Waxman 2001; Panaitescu & Kumar 2001; Berger et al. 
2003a). The price of accepting this picture would be to 
abandon the current afterglow theory completely, which 
has been well tested and proven adequate to deal with the 
bulk of the afterglow data. If, on the other hand, we take 
£■52 ~ 1 as inferred from data (e.g. Berger et al. 2003a), 
this narrow-beam picture would imply that GRBs have a 
total energy of order 10^^ ergs. It would also imply that 
the number of GRBs is similar to the number of SN Ib/c, 
whereas a sample of SN Ib/c reveals only 3% have radio 
afterglows (Berger et al. 2003c). 

A variant of the uniform jet model is to interpret GRBs 
as on-beam detections and XRFs as off-beam detections 



(Yamazaki, loka & Nakamura 2003a). This model may 
have the prospect of both interpreting the correct XRF-to- 
GRB ratio and preserving the standard afterglow model, 
but still requires GRB jets to have a large dispersion of 
opening angles. Furthermore, the XRFs' afterglows in 
such a model should resemble those in the "orphan af- 
terglows", i.e., initially rising and peaking at a time tpk 
when 1/r = 9y, where F is the Lorentz factor within the 
uniform jet cone, and 9v is the viewing angle (e.g. Granot 
et al. 2002). Since 9y > 9j in this model, one should ex- 
pect that tpk is typically larger than the typical tj in GRB 
afterglows. The recent afterglow observations for XRF 
030723 indeed show an initial rising lightcurve (Huang et 
al 2003, and references therein), but tpk is around 0.1 day, 
much smaller than the typical tj for GRBs, which is typ- 
ically several days. In order to interpret the XRF as an 
off-beam GRB, the jet opening angle has to be anoma- 
lously small. The lightcurve, however, is consistent with 
the standard on-beam afterglow with the peak being due 
to crossing of the typical synchrotron frequency across the 
optical band (e.g. Kobayashi & Zhang 2003). Similarly, 
the optical afterglow data for XRF 020903 (Soderberg et 
al. 2003) indicate that the lightcurve already decays start- 
ing from 0.9 days, which is not consistent with an orphan 
afterglow lightcurve. 

3. A QUASI-UNIVERSAL GAUSSIAN-LIKE STRUCTURED 
JET MODEL 

Below we will argue that the current GRB-XRF prompt 
emission and afterglow data are compatible with a model 
in which the GRB-XRF jets have a riuasi-univcrsal struc- 
ture, where the jet energy distribution is axially symmet- 
ric and Gaussian-like, i.e., a jet which still has a typical 
angle, with a mild variation of energy within this jet typ- 
ical angle and a rapid (e.g. exponential) decrease of en- 
ergy outside the typical angle (Zhang & Meszaros 2002a; 
Lloyd- Ronning et al. 2003). For analytical purposes, we 
approximate the angular distribution of the jet energy as 



e(^) = eoe 



(3) 



The initial Lorentz factor should also have an angular de- 
pendence, and we take it as a free function, since it is in 
principle independent and since it does not directly infiu- 
ence the GRB-XRF population study ^. The motivation 
for introducing a Gaussian-like jet are dual. First, it pre- 
serves a characteristic angle for the jet, which is more 
consistent with the jet structure generated in the numer- 
ical simulations for the coUapsar model (Zhang, Woosley 
& MacFadyen 2003a; Zhang, Woosley & Heger 2003bj. 
It also alleviates the divergence problem of a power-law 
model at small angles by introducing a smooth functional 
profile within the characteristic jet angle. Second, at larger 
angles (beyond 9o), the energy decrease is steeper than the 
power law model, so that in order to get an XRF whose 
Eiso is (10^ — 10^) times lower than the typical GRB Eiso, 
one only needs to have a viewing angle 9^ {3 — 4)^o- 
This greatly reduces the predicted number of XRFs. 



®The discussions here are relevant when the prompt emission energy is roughly proportional to the kinetic energy, which is verified at least 
for GRBs since both geometry-corrected energy components are roughly standard. 

'^It, however, influences the rising time of the afterglow lightcurve. See below for further discussion. 
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We perform a simple Monte Carlo simulation to ver- 
ify this ansatz. Since it is unreasonable to expect that 
all GRB-XRFs are exactly the same, we allow some scat- 
ter of the model parameters. Such scatter is also needed 
(Lloyd-Ronning et al. 2003) to reproduce the observed 
Eiso — dj correlations from the afterglow data (Frail et al. 
2001; Bloom et al. 2003), even for the k = 2 power-law jet 
model. We randomly generate 10000 bursts, each of which 
has a jet structure in the form of eq.(3). We define the to- 
tal energy within the jet as Ej = 2tt JJ^^^ e(^) sin^rf^. For 
small angles, this is 27reo^?o (Zhang & Meszaros 2002a), 
but here we use a more rigorous integration which is also 
applicable for large angles. We assume that the distri- 
butions of both Ej and Oo are lognormal. For each re- 
alization, we use the generated Ej and to derive eo 
and to obtain the jet structure according to (eq.[3]). For 
each burst, we also generate a viewing angle with the 
probability P{6y)d9y cx sm{6y)d6y. The corresponding 
4:Tre{9v) is then assigned to be the E^go for that partic- 
ular burst. The cosmic rest frame Ep of that burst is 
generated via eq.(l) with a lognormal variation (with a 
standard deviation of 0.3), and the observed peak energy 
is Ep{ohs) = Ep/{1 + z). The redshift distribution of the 
bursts is assumed to trace the cosmic star-forming rate 
(Rowan- Robinson 1999). Standard cosmological parame- 
ters arc adopted, i.e., Hq = 70 km s^^ Mpc^^, n„i = 0.3 
and l^A = 0.7. No redshift-evolution for the burst pa- 
rameters are assumed. We then assign a redshift for each 
burst, and calculate the distance and the energy fluence 
for that burst. Finally we place a fluence truncation of 
5 X 10^^ erg cm^^ to reflect the instrument sensitivity 
limit. This number matches the faintest HETE-2 XRF 
(Fig.l of Lamb et al. 2003), and we regard it as reflecting 
HETE-2 sensitivity. (We have experimented with chang- 
ing the fluence truncation limit. In general, the total num- 
ber of the detectable GRB-XRFs increases with a deeper 
fluence truncation, and the fraction of XRFs from the 
whole population increases mildly.) The simulated bursts 
are plotted in the £'p(o6s)-fluencc plane (Fig.l), together 
with the BeppoSAX and HETE-2 prompt emission data. 

The model can be also compared to the afterglow light 
curve break data (Bloom et al. 2003). In a structured jet, 
if the jet structure is steep enough (e.g. for the power law 
jets), the viewing angle defines the jet break time (Rossi 
et al. 2002; Zhang & Meszaros 2002a; Kumar & Granot 
2003; Granot & Kumar 2003; Panaitescu & Kumar 2003; 
Wei & Jin 2003; Salmonson 2003). For a Gaussian jet, 
although the viewing angle defines the jet break time for 
Oy > ^Oj the jet structure is only mild within the typical 
angle ^Oj so that for 9y < Oq, it is rather than 9y which 
defines the jet break time (Kumar & Granot 2003). In our 
simulation, we assign a jet break angle 9j = max(9o,9-u) 
for each burst^. The fluence-truncated bursts are plotted 
in the Eiso — 0j plane (Fig.2), with the afterglow data. 

Figures 1 and 2 show our simulation results. The pa- 
rameters adopted are \og{Ej/l erg) = 51.1 ± 0.3 and 
log(0o/l rad) = — 1.0±0.2, where the error is the standard 
deviation of the lognormal distribution. This corresponds 
to a quasi-standard explosion energy of Ej = 1.3lJ;g x 10^^ 
ergs, and a quasi-standard jet configuration with the typ- 



ical opening angle of 6*0 = '^■'^-21 degrees. This is in 
excellent agreement with coUapsar simulations (Zhang et 
al. 2003a). Figure 1 shows that this model is com- 
patible with the prompt emission data collected by Bep- 
poSAX and HETE-2. Defining XRFs as those events with 
Ep{obs) < 25 keV (Soderberg et al. 2003), the simulated 
number of XRFs is roughly 1/3 of the whole GRB-XRF 
population within the fiuence-truncated observed sample, 
as in the observations. Figure 2 shows that the same set 
of parameters is compatible with the jet angle data within 
the framework of the standard afterglow model. 

Among 10000 simulated bursts, 702 GRB-XRFs survive 
the fluence-selection effect. Since this takes viewing angle 
effects into account, the simulation suggests that the true- 
to-observed number ratio for the entire combined GRB- 
XRF population is about 14 (or 21 for GRB plus X-ray- 
rich bursts, or 42 for GRBs alone). The ratio should 
decrease for higher sensitivities. This number is much 
smaller than the beaming correction factor (~ 500) in the 
uniform jet model for the GRB population (Frail et al. 
2001; van Putten & Regimbau 2001), suggesting that the 
number of GRB-XRF progenitors required is smaller than 
previously thought. 

Kumar & Granot (2003) have performed a numerical 
hydrodynamical modeling of the evolution of a Gaussian- 
jet, which reveals a substantial energy redistribution in 
the jet structure during its evolution. The resultant 
lightcurvcs (their Fig. 5a) are consistent with most of the 
jet break data. For 9y < ^^o, the lightcurvcs are simi- 
lar to those of a uniform jet model which is consistent 
with the data. For 9v > 9o, & jet break is visible around 
the time of r(6'^,) ^ 9~^. In their calculations, the struc- 
ture of T{9-u) is also taken as a Gaussian, so that when 
9v ^ 9o, the lightcurves become analogous to those of or- 
phan afterglows characterized by a late rising. However, 
if the initial Lorentz factor at large viewing angles is high, 
i.e., ro(^^t,) 9~^, the rising segment of the lightcurve 
would be shifted to earlier times, followed by a decaying 
lightcurve (Zhang & Meszaros 2002a). For 9^ » ^o, a 
lightcurve bump may be expected for 9^ » ^0 when the 
jet axis becomes visible (J. Granot, 2003, personal commu- 
nication), since the main power at the jet axis is initially 
not in the relativistic light cone. However, the energy re- 
distribution effect (e.g. Fig.l of Kumar & Granot 2003) 
may smear the bump to be less significant. In any case, 
XRFs in the current model only imply 9y to be at most 
4:9o- Observationally, light curve breaks are currently well 
monitored mainly for those bursts whose break time is 
typically days, which corresponds to small jet angles (i.e. 
9y<9o). The afterglow data for the recent XRF 030723 
(Huang et al. 2003 and references therein), on the other 
hand, shows a lightcurve re-brightening after 10 days. For 
an initial Gaussian jet (in the prompt phase), the jet struc- 
ture would evolve during the afterglow phase because of a 
pole-to-equator energy flow. For a large viewing angle cor- 
responding to an XRF, the late-time afterglow energy in 
the viewing direction could be much larger than the energy 
in the prompt phase. This is consistent with the recent ra- 
dio afterglow observation for XRF 020903 (Soderberg et 
al. 2003). We therefore conclude that the model discussed 



*Thc result docs not differ too much from the present simulation if we assumed that 9,j is exclusively the agent which defines the jet breaJf 
(e.g. Lloyd-Ronning et al. 2003). The reason is that the probability for small viewing angles (which maJse the difference) is small. 
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here is not inconsistent with ah the current GRB-XRF 
afterglow data. 

4. DISCUSSION 

Wc have argued that in order to incorporate both the 
prompt emission and afterglow data for GRBs and XRFs 
within a unified theoretical framework, both the simple 
fc = 2 power-law universal jet model and the on-bcam uni- 
form jet model encounter difficulties. With a Monte Carlo 
simulation, we show that the current data are compatible 
with a quasi-universal Gaussian-like structured jet model. 

Other possible models may still interpret the data. 
First, a uniform jet also produces an equivalent isotropic 
energy outside the jet cone due to the off-beam Doppler 
effect. So an off-beam model for XRFs (Yamazaki et al. 
2003a) does not over-generate XRFs. However, the de- 
cline of the isotropic energy at larger viewing angles is 
even steeper than exponential (e.g. Yamazaki, Yonetoku 
& Nakamura 2003b). So one can only generate XRFs at 
viewing angles slightly larger than the jet opening an- 
gle. This tends to under-generate XRFs. The afterglow 
issue as discussed in §2 also poses some constraints on 
this model. In any case, a population simulation like the 
one presented here is needed to validate that model. Sec- 
ond, it is possible that some or even all GRB-XRF jets 
include two (or more) components (Li & Chevalier 1999; 
Frail et al. 2000; Lipunov, Postnov & Prokhorov 2001; 
Berger et al. 2003b; Sheth et al. 2003; Huang et al. 



2003), as reproduced from the coUapsar models (Zhang et 
al. 2003a, 2003b; Mcszaros & Recs 200f; Ramirez-Ruiz, 
Celotti & Rees 2003). In such a model, the jet structure 
may be approximated as a superposition of two (or more) 
uniform or Gaussian components. Such a model may be 
also able to reproduce the GRB-XRF data. However, the 
afterglow lightcurves could show distinct signatures un- 
der certain conditions (Berger et al. 2003b; Huang et 
al. 2003). The lack of such signatures in most of the 
afterglow lightcurves poses constraints on the parameters 
for such two-component structured jet models. Third, 
we do not exclude the possibility that ;i k = 2 power 
law structure extends only up to a certain angle, above 
which the jet energy has a much steeper decline (Zhang & 
Meszaros 2002a). However, the current Gaussian model 
is the simplest among these, and can accommodate the 
widest range of prompt emission and afterglow data for 
the majority of GRBs and XRFs. More detailed Monte 
Carlo simulations for this as well as other models and their 
comparison to a wider set of detailed prompt emission and 
afterglow data will be presented elsewhere (Dai & Zhang 
2003, in preparation, see also Liang, Wu & Dai 2003). 
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01128, and the MoneU Foundation. 



REFERENCES 



Amati, L. et al, 2002, A&A, 390, 81 

Berger, E., Kulkarni, S. R. & Frail, D. A. 2003a, ApJ, 590, 379 
Berger, E. et al. 2003b, Nature, 426, 154 

Berger, E., Kulkarni, S. R., Frail, D. A. & Sodcrbcrg, A. M. 2003c, 
ApJ, 599, 408 

Bloom, J. S., Frail, D. A. & Kulkarni, S. R. 2003, ApJ, 594, 674 
Frail, D. A. ct al. 2000, ApJ, 538, L129 
. 2001, ApJ, 562, L155 

Frecdman, D. L. & Waxman, E. 2001, ApJ, 547, 922 
Granot, J. & Kumar, P. 2003, ApJ, 591, 1086 

Granot, J., Panaitescu, A., Kumar, P. & Woosley, S. E. 2002, ApJ, 
570, L61 

Harrison, F. A. et al. 1999, ApJ, 523, L121 
Heise, J. et al. 2003, AIP Conf.Proc. 662, 229 

Huang, Y. F., Wu, X. F., Dai, Z. G., Ma, H. T. & Lu, T. 2003, ApJ, 

submitted (astro-pli/0309360) 
Kippen, R. M. et al. 2003, AIP Conf.Proc. 662, 244 
Kobayashi, S. & Zhang, B. 2003, ApJ, 582, L75 
Kulkarni, S. R. et al. 1999, Nature, 398, 389 
Kumar, P. & Granot, J. 2003, ApJ, 591, 1075 

Lamb, D. Q., Donaghy, T. Q. & Graziani, C. 2003, to appear in proc. 
2nd VERITAS Symposium on TeV Astrophysics, Chicago, Illinois 

(astro-ph/0309456) 
Li, Z.-Y. & Chevalier, R. 1999, ApJ, 526, 716 
Liang, E. W., Wu, X. F. & Dai, Z. G. 2003, MNRAS, submitted 
Lipunov, V. M., Postnov, K. A. & Prokhorov, M. E. 2001, Astron. 

Rep. 45, 236 



Lloyd, N. M., Petrosian, V. & Mallozzi, R. S. 2000, ApJ, 534, 227 
Lloyd-Ronning, N. M., Dai, X. & Zhang, B. 2003, ApJ, in press 

(astro-ph/0310431) 
Mcszaros , P. & Rees, M. J. 1997. ApJ, 482, L29 
. 2001, ApJ, 556, L37 

Panaitescu, A. & Kumar, P. 2001, ApJ, 560, L49 
. 2003, ApJ, 592, 390 

Ramirez-Ruiz, E., Celotti, A. & Recs, M. J. 2002, MNRAS, 337, 1349 
Rhoads, J. E. 1999, ApJ, 525, 737 

Rossi, E., Lazzati, D. & Rees, M. J. 2002, MNRAS, 332, 945 
Rowan-Robinson, M. 1999, Ap&SS, 266, 291 
Sakamoto, T. et al. 2003, ApJ, submitted (astro-ph/0309455) 
Salmonson, J. 2003, ApJ, 592, 1002 

Sodcrbcrg, A. M., Kulkarni, S. R., Berger et al. 2003, ApJ, submitted 

(astro-ph/0311050) 
Sheth, K. ct al. 2003, ApJ, 595, L33 

van Putten, M. H. P. M. & Regimbau, T. 2003, ApJ, 593, L15 
Wei, D. M. & Jin, Z. P. 2003, A&A, 400, 415 

Woosley, S. E., Eastman, R. G. & Schmidt, B. P. 1999, ApJ, 516, 
788 

Yamazaki, R., loka, K. & Nakamura, T. 2003a, ApJ, 593, 941 
Yamazaki, R., Yonetoku, D. & Nakamura, T. 2003b, ApJ, 594, L79 
Zhang, B. & Mcszaros , P. 2002a, ApJ, 571, 876 
. 2002b, ApJ, 581, 1236 

Zhang, W., Woosley, S. E. & MacFadyen, A. I. 2003a, ApJ, 586, 356 
Zhang, W., Woosley, S. E. & Heger, A. 2003b, ApJ, submitted (astro- 
ph/0308389) 



5 




Fig. 1. — The -Ep(o6s)-fluence diagram of the simulated GRB-XRFs (small dots) as compared with the HETE-2 and BeppoSAX data 
(heavy dots with error bars, data from Lamb et al. 2003). 
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Fig. 2. — The Ejso — Sj diagram of the simulated GRB-XRFs (small dots) as compared with the afterglow jet break data (heavy dots with 
error bars, data from Bloom et al. 2003). 



